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ABSTRACT 

Linear instability of two equilibrium configurations with either poloidal (I) or toroidal 
(II) dominant magnetic field components are studied in thin vertically-isothermal Ke- 
plerian discs. Solutions of the stability problem are found explicitly by asymptotic 
expansions in the small aspect ratio of the disc. In both equilibrium configurations 
the perturbations are decoupled into in-plane and vertical modes. For equilibria of 
type I those two modes are the Alfven-Coriolis and sound waves, while for equilibria 
of type II they are the inertia-Coriolis and magnetosonic waves. Exact expressions for 
the growth rates as well as the number of unstable modes for type I equilibria are 
derived. Those are the discrete counterpart of the continuous infinite homogeneous 
cylinder magnetorotational (MRI) spectrum. It is further shown that the axisymmet- 
ric MRI is completely suppressed by dominant toroidal magnetic fields (i.e. equilibria 
of type II). This renders the system prone to either non-axisymmetric MRI or non- 
modal algebraic growth mechanisms. The algebraic growth mechanism investigated in 
the present study occurs exclusively due to the rotation shear, generates the inertia- 
Coriolis driven magnetosonic modes due to non-resonant or resonant coupling that 
induces, respectively, linear or quadratic temporal growth of the perturbations. 
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1 INTRODUCTION 

The stability properties of Keplerian disks have been a focus of intensive investigation of theoretical astrophysicists over the last 
decades, pertaining to the problem of angular momentum transfer in accretion disks. Traditionally, two routes have been taken 
by researchers in order to elucidate various aspe cts of that topic . Thus, the results of the analytical study of magneto-rotational 
instability (MRI) in an infinitely long cylinder Velichov 1 19591 ): Chandrasechar ] |l960l)l have b e en ad opted for the thin disks 
in order to derive criteria for the spectral stability under various conditions [Balbus fc Hawlev ( 199 ll )]. On the second route, 
numerical calculations have been performed, usually by employing the shearing-box model, in or der to study the nonlinea r 



evolution of small pe r turbations, and th e ir subsequent development into full sustainable turbulence Brandenburg et alj (| 19951 ); 
Hawlev et all |l99rj ): lBodo et al.l |200sl ): lRegev fc Umurhanl |200sl )]. The main thrust of both routes has been to establish the 



MRI as the main generator of sustainable turbulence. Notwithstanding accumulating analytical as well as numerical results, 
the notion of the MRI as a driver of turbulence and angular momentum tr ansport has r ecent l y been shown to suff e r from 
non-trivial difficulties in the sheer use of the shearing box in the simulat ions Bodo et al.l J200ct); Regev fc Umurhan I (120081)1 
as wel l as in doubts regarding the numerica l conve rgence and resolution Lesur fc Longarettil (|2007h ; Fromang fc Papaloizo 



I 2007"); iFromang et al.l (|2007l ); |Pessah et al.l (|2007h ]. In addition, it has recently been shown that if the thin disk geometry is 
taken into account both the grow th rates as well as the number of unstable MRI m odes are greatly reduced and are decreasing 
functions of the disk thickness [ Coppi fc Keves 1 20031 ); Liverts fc Mon d (2009)]. In parallel, an alternative route that has 
recently emerged within which the dynamical response of the thin disks is investigated asymptotically in the small thickness- 
to-radius ratio has been shown to lead to fruitful and physically sound results. In particular, applying such strategy it has 
been shown that algebraic non-modal instabilities are the primary source of pure hydrodynamical activity Umurhan et al 
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l 2006tl : lRebusco et all (|2009l ); lshtemler et all j2O10h ]. Furthermore, the mechanism for such non-exponential g rowth has been 



Shtemler et al 



identi fied as the resonant as well as non-resonant interaction between inertia-Coriolis modes and sound waves 
1 201(3)]. Interestingly enough, such an investigation has never been carried out for magnetized disks that are described by the 
magnetohydrodynamic (MHD) model. It is thus one of the main goals of the current work to carry out a comprehensive study 
of the MHD response of thin Keplerian disks taking into account their true thin-disk geometry. 

In that connection note that although the magnetic field configuration in real discs is largely unknown, in the initial 
stage of the disc rotation the poloidal field is commonly accepted as the dominant component. However, observations and 
numerical simulations indicate t hat toroidal magnetic fie l d may be of the same order of magni t ude, or even dominate th e 



poloida l mag netic field [see e.g. iTerquem fe Papaloizou I (| 199 



6): 



Proga I (20Cy )]. In particular, the unstable modes give rise to MRI 



Papaloizou fc Terqueml Jl997h ; lHawlev fc Krolik I (|2002h : 



Balbus fc Hawlev ( 199 ll )] which amplifies considerably th e 



toroidal component o f the magnetic field to the point that it may dominate the original axial field Pessah fc Psaltia ( 2005 ); 



Begelman fc Pringld ([20071 ')] . It is of great importance therefore to investigate how the magnetic field topology influences the 



stability properties of the disc. The way to do that is to follow the linear problem in terms of a toroidal dominated magnetic 
field instead of poloidal one leading to MRI. Thus, one more goal of the present work is investigation of two types of magnetic 
field equilibrium configurations with either poloidal or toroidal dominated magnetic field component and their relations with 
spectral and non-modal mechanisms of instability in thin discs. 

The paper is organized as follows. The dimensionless governing equations and their approximation to leading order in 
small aspect ratio of the steady-state disc are presented in the next Section. Two magnetic field configurations with dominated 
axial (I) and dominated toroidal (II) equilibrium magnetic fields are also described in that Section. Sections 3 and 4 describe 
the principal hydro-magnetic modes which can propagate in thin Keplerian discs for I and II types of equilibria, respectively. 
Summary and discussion are presented in Section 5. In Appendix A the general relations for comparable in-plane components 
of the equilibrium magnetic field are discussed in the limit of large plasma beta, and the qualitative stability analysis is made. 



2 THE PHYSICAL MODEL FOR THIN KEPLERIAN DISCS 

The stability of radially as well as axially stratified rotating plasma in thin vertically isothermal discs threaded by a magnetic 
field is considered. Viscosity, electrical resistivity, and radiation effects are ignored. 



2.1 Governing equations 

As a first step, all phy si cal variables are tran sformed to non-dimensional quantities by using the following characteristic values 



Shtemler et all |2009h ; Ishtemler et al.l (|20ich ] 



1 ?"* 
t* = — , V» = — , L* — V*t t , m* = m,i, n* 



= V, 2 , cs« = VT*/m«, P. = m»n«c s 2 , j« = *-, E* = -=-=.. (1) 

Here fi* = (GMc/r*) 1 ' 2 is the Keplerian angular velocity of the fluid at the characteristic radius r* that belongs to the 
Keplerian portion of the disc; G is the gravitational constant; M c is the total mass of the central object; c is the speed of 
light; $* is the characteristic value of the gravitational potential; the characteristic mass and number density equal to the ion 
mass and number density, m» = rrn and n* = rij. The characteristic values of the electric current density and electric field, 
j» and E, have been chosen consistently with Maxwell's equations; cs* is the characteristic sound velocity; T» = T(r*) is the 
characteristic temperature; the characteristic magnetic field is specified below depending on the magnetic field configuration: 
P* = B z (r t ) if the axial equilibrium magnetic field dominates or of the order toroidal one, and B, — Bg(r t ) otherwise; 
the dimensional equilibrium temperature and the two components of the equilibrium magnetic field, T(r) and B z (r), Bg(r), 
respectively, are free functions. 

The resulting dimensionless dynamical equations for vertically isothermal discs are: 

DV 1 VP „^ 1 j x B 
— V<1> -I 

Dt MJ n /9M| 

dn 

In 
<m 
~dt 

E = V x B, (5) 



79 ^® + ^T79- , ( 2 ) 

+ V ■ (nV) = 0, (3) 
V x E = 0, V ■ B = 0, (4) 



P = nT. (6) 
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Here VP = c|Vn for vertically isothermal discs, and the dimensionless equilibrium sound speed is given by c| = dP/dn = 
T(r). Standard cylindrical coordinates {r, 8, z} are adopted throughout the paper with the associated unit vectors {i r , ie, i z }; 
V is the plasma velocity; t is time; D/Dt = d/dt + (V ■ V) is the material derivative; <E>(r, z) = — (r 2 + z 2 )' 1 ^ 2 is the 
gravitational potential due to the central object; B is the magnetic field, j = V x B is the current density; E is the electric 
field; P = P e + Pi is the total plasma pressure; Pi — niTi are the partial species pressures (I = e, i); T = T e = T» is the plasma 
temperature; subscripts e and i denote electrons and ions, respectively. Note that a preferred direction is tacitly defined here, 
namely, the positive direction of the z axis is chosen according to positive Keplerian rotation. The dimensionless coefficients 
Ms and /3 are the Mach number and the characteristic plasma beta, respectively: 

cs* ' P B 2 

Zero conditions at infinity for the in-plane magnetic field are adopted, namely: 



M S = -^, /3 = -l. (7) 



B r = 0, B = for z = ±oo. (8) 

The boundary condition for the axial velocity will be specified below as the least possible divergence at infinity. As will be 
seen later on, despite such unbounded growth the axial mass flux tends to zero far away from the mid-plane due to the fast 
vanishing of number density at z — ±oo. To simplify the further treatment of Maxwell's equations, both the hydro-magnetic 
basic configuration and perturbations are assumed to be axisymmetric. 



A common property of thin Keplerian discs is their highly compressible motion with large Mach numbers [ iFrank et al 



1 20021 )] . Furthermore, the characteristic effective semi-thickness of the equilibrium disc H„ — H(r») (H — H(r) is the local 
semi-thickness) is defined so that the disc aspect ratio e equals the inverse Mach number: 

Thus, the thin disc approximation means 



^ = V§t= e<<1 ' (10) 

where 



Ms = ^, K=rA, n.=J^, c s * = Vn, H* = ^ . 



cs* V ri ll* 

The smallness of e means that dimensionless axial coordinate is also small, i.e. z/r* ~ e (\z\^H,), and consequently the 
following rescaled quantities may be introduced in order to further apply the asymptotic expansions in e [similar to Shtemler 
et al. (2009); Shtemler et al. (2010)]: 

C=f~ £ °, ff(r) = ^~ e °, (11) 
where H(r) = cs(r)/Cl(r) is the scaled semi-thickness of the disc. 

2.2 Equilibrium configurations 

We start by deriving the steady state solution. It is first noted that the asymptotic expansion for the time-independent 
gravitational potential is given by: 

S(r,C) = Hr) + e 2 ^(r,0, *(r) = --, 4>{r, fl = k^ 2 (r) + 0(e 2 ), r > 1 » e. (12) 

r 2 

Substituting (|12|l into Eqs. (IS}-© and setting the partial derivatives with respect to time to zero yield to leading order in e 

Yl = dmr± c 2 s (r)dn_ d$(r,Q 
r dr ' n d( d( ' 1 ' 

Thus, the velocity as well as the number density are as follows: 

K = o(e 2 ), Vg = e Ve{r) + O(e 2 ), V z = o(e 2 ), n S e°n = e° N(r)v{r,), (14) 
where o(e) -C e, 0[e) ~ e, and 

%(r) = rfl(r), n(r) = r~ 3/2 , P(r,) = exp(-r? 2 /2), n = t/H(r). (15) 

Two quite different equilibrium magnetic configurations will be considered below, namely one with comparable magnitudes 
of the axial and toroidal components of the magnetic field (mainly studied for the dominant poloidal magnetic field, while 
the general case of comparable components of the magnetic field is discussed in Appendix A), and the other with dominant 
toroidal component. The two equilibria (denoted I and II below) are distinguished by different scaling of the physical variables 
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Table 1. Scales in e of the equilibrium variables for two types of equilibria (I) and (II). 



/ = e S f 


n 


V r 


V e 


v z 


13, ■ 


Be 


B z 


jr 


39 


jz 


I, S 





> 2 





> 2 


> 2 








> 2 








II, s 





> 2 





> 2 


> 2 





1 


> 2 


1 






with e which generally may be written as /(r, £) = e S f(r,C)- All equilibrium variables are written in the leading order in e, 
and depend on the radial variable only. The exceptions are the number density and the pressure that depend on the axial 
coordinate in a self-similar manner with radius-dependent amplitude. The toroidal and axial magnetic fields as well as the 
disc thickness and the amplitude factor, N(r), in the number density are arbitrary functions of the radial variable. Those 
functions specify the equilibrium state. 

To start the equilibrium description it is first assumed that both types of equilibria under the current study are char- 
acterized by a toroidal component of the magnetic field that is of order e°. In contrast, in one equilibrium system (I) the 
axial component of the magnetic field is of order e° as well, while in the other system (II) it is of order e. These assumptions 
together with relations (|12p - (|15p . determine the order in e of the rest of the physical variables: 

(I) Comparable equilibrium magnetic field components 

B r = o(e 2 ), Bg^e°B (r), B z = e° B z (r), j r = o(e 2 ), jg = e°j e (r) 

(II) Dominant equilibrium toroidal magnetic field 

B r = o{e 2 ), B e ^e°B e (r), B z St eB z (r),j r = o(e 2 ), j e ~ eje(r) = 

For convenience the results are summarized in Table 1. It is noted finally that to lowest order in e the magnetic field 
configurations under consideration do not influence the steady-state properties of the disk. As will be seen in the following 
sections, this situation changes dramatically when small perturbations are considered. 



dB z . o- oldirBe) 
dr r dr 



(16) 



dB z . _ ld{rBg) 

-e—r- , Jz = e = e ^ • 

dr r dr 



(17) 



2.3 Perturbed thin discs 

In general for the unsteady nonlinear case the dependent variables are scaled in e in the following way: 

f(r,C,t) = e s f(r,0+e s 'f'(r,C,t). (18) 

Here f(r, C, t) stands for any dependent variable, the bar and the prime denote equilibrium and perturbed variables; each 
perturbed variable is characterized by some power S' that is different for two types of equilibria I and II. The various cases 
are summarized in (Table 2). 



3 DYNAMICAL EQUATIONS FOR THE PERTURBED DISC: TYPE I EQUILIBRIA 

Consider the general dynamic equations for the perturbed disc equilibria of type I, namely thin discs under the effect of 
a magnetic field with comparable poloidal and toroidal components. Then the detailed stability study is carried out for 
zero toroidal magnetic field, while a qualitative analysis of the comparable poloidal and toroidal components is presented in 
Appendix A. 



3.1 The reduced nonlinear equations 

Substituting (|12[) - (|18[ ) into J2| - ([6]) yields to leading order in e the following system of the reduced non-linear equations that 
governs the dynamical behavior of thin discs: 

^-20(^ + ^-1^1^=0, (19) 
ot oQ p n + n' oQ 

ev[ + i d£fi) v , a^iB^e*, 

Ot r dr oC, p n + n' 

dV' z _ 2 n 8 n' ,l 6M 8, _g_ > 1 8 /2 1 B' B 2 + B r \ 
ot n + n' o£ n p oQ n + n' 2o( p n + n' 

On' d[{n + n')V:] _ 

8t + d( ( > 



Spectral and algebraic instabilities in Keplerian discs 

Table 2. Scales in e of the perturbed variables for two types of equilibria (I) and (II). 



f = e Sf + e S f 
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V r 


Vg 


v z 


Br 


Be 


B z 


jr 39 jz 


I, S' 





1 


1 


1 








1 


-1 -1 


II, S' 











1 








1 


-1 -1 



The system of equations ()19[) - (|24|) is subject to the boundary condition of least possible divergence for the axial velocity at 
infinity along with the conditions for the in-plane magnetic field as follows from (8) 

B' r = 0, B' e = for C = ±°o. (25) 

Both poloidal and toroidal components of the perturbed magnetic field, B' r and B' z , are expressed through the magnetic flux 
function, ty' , which provides the divergent free magnetic field, V B = 0. Here the equation for B' r is obtained by differentiating 
by £ the equation for ^P', while the corresponding equation for the perturbed axial magnetic field, B' z , decouples from the rest 
equations, and drops out from the governing system of equations (|19|l - (|24|l . Equation (|21[) has been derived with the help of 
the second steady-state relation (|14|) . Furthermore, the radial derivatives drop out from the system of Eqs. (|19[) - ()25|) . This, 
it should be stressed, is so not due to a frozen coefficients assumption, but is a direct result of the thin disc geometry. An 
important conclusion from the absence of radial derivatives in the reduced system of equations of the thin disc approximation 
is that transient non-exponential growth is not possible in such configuration. As will subsequently be shown (see Section 4 of 
the present paper), this picture changes dramatically when the dominant magnetic fields are toroidal, in which case the only 
possible axisymmetric instability is due to non-modal algebraic transient growth. 

As the radial coordinate is a mere parameter in the set of the reduced equations, it is convenient to replace the physical 
variables by the following self-similar quantities: 

r = t, r,= ^- y (26) 

such that the derivatives in the new and old variables are related as follows: 
1 



Ot Or' OC H(r)&n 

Finally, system (|19[) - (|25[) may be written in a simpler form by introducing the following scaled dependent variables: 

V n' B' 

v(T,r,v) = — 77, v(r,r,r,) = b(r, r,rj) = -g-yr-- (28) 

cs(r) N(r) B z (r) 

Below for convenience and with no confusion the notation t for the time variable is reinstated instead of the new variable r. 
This yields the following system of equations that depend parametrically on the radius: 

1 dv r 1 1 0b r _ OVr . , 

fl(r) dt Pz(r)Hv) + v dr, ~ Vz Or, ' ( > 

1 Ovg 1 1 1 0b e 0v e , on , 

+ -V r — - — — = —V z —— , (30) 



Q(r) Ot 2 Pz{r)v{r)) + v Or, drj 
1 0v z+ Hv) d_ ( ^ ) + £^l± [ be ] = _l± [v 2 + J: j + (31) 



O(r) Ot v{rj) + vdr) v(rf) /3 z (r) Or/ v{rj) + v 2 Or, f3 z (r) v{r,) + v 

1 dv + d[{v{Tj)+v]v z =Q| 



Q{r) Ot Or, 

1 0b r 0v r _ 0(v z b r ) 

O(r) Ot Or, Or, 



(33) 



1 Obg Ovg 3 0{v z bg) =. 0v z 

-K~ + n b r = HZ S ( r )^T' ( 34 ) 



O(r) Ot Or, 2 Or, v ' Or, 

supplemented by the boundary condition of least possible divergence for the axial velocity along with the vanishing conditions 
for the in-plane magnetic field at infinity. Here the Keplerian epicyclical frequency, x( r ) = ^( r ) nas been employed; v(rj) is 
the scaled equilibrium density; j3 z (r) and fie(r) are the local axial and toroidal plasma beta functions: 
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BI(r) B^(r) y /3 e (r) B z (r) 

where both local parameters /3 z (r) and Pe(r) are proportional to the characteristic plasma beta, /3. 

Relations (|29[) - (|35|) form the full nonlinear MHD problem in the thin disc approximation and are named as defined above, 
the reduced nonlinear equations. 

3.2 The linear problem. 

Assuming now that the perturbations are small, the system of equations (|29[) - (|34|1 may be linearized about the steady-state 
equilibrium solution. This yields: 

1 dv r 1 db r 

- - — = 0, (36) 



Q{r) dt Mr)Hv) 9v 

1 dv g 1 1 db e , 

+ ■xVr - - -jr- = 0, (37) 



0(r) dt 2 &(r)P(r?) d V 
1 db r dv. 



fi(r) dt dr) 



= 0, (38) 



1 dbe dv e + K = _ S{r) dv lj (39) 



O(r) at dr) 2 y J drj 

1 9 ^ + JL(J^) = -pp-JL[-h-], (40) 



O(r) 9t 1/(77) <3z( r ) a?? ^(t?) 

1 a, + a[£^K] =a 



n(r) at a?? 

For definiteness the characteristic magnetic field is specified as the dimensional axial component of the equilibrium magnetic 
field B, = B*(r„). 

The general case of comparable poloidal and toroidal components of the equilibrium magnetic field is presented in 
Appendix A for perturbations of type I. A qualitative analysis is done of the influence of the equilibrium toroidal magnetic 
field on the disc stability in the limit of large plasma beta. In the case of small plasma beta (/3 < 1) two principal modes, 
namely the Alfven-Coriolis (AC) and the magnetosonic (MS) modes, are strongly coupled, while in the limit of large plasma 
beta, AC mode decouples on its characteristic scale from the MS mode. This leaves the resulting dispersion relation the same 
as in the case of zero toroidal magnetic field, and the influence of the toroidal magnetic field is reduced to excitation of the 
AC-driven MS mode. On the characteristic scale of MS's mode, a stable MS mode decouples from the AC mode and the 
influence of the equilibrium toroidal magnetic field is reduced to excitation of the stable MS-driven AC mode. 

To simplify the analysis, the equilibrium toroidal magnetic field is set to zero for type I equilibria, or equivalently, 
/3e(r) = 00. Under such conditions the reduced linearized system (|36[) - (|41|l is divided into two decoupled sub-systems that 
describe the dynamics of two different modes, namely: the Alfven-Coriolis and the sound modes. It is indeed the approximation 
that is adopted in the remaining of the current analysis of type I equilibria. 

3.3 Linear stability analysis for the Alfven-Coriolis modes. 

We start by representing the perturbations up to a radius-dependent amplitude factor as follows: 

f(r, 77, t) = exp[-iX(r)t]f(r, r,), (42) 
where A(r) is the complex eigenvalue 

A = A + tT. (43) 

Substituting (|42I) - (I43I ) into H36I)- (I41I ) results in the following system of linear ordinary differential equations for the 
perturbed velocity as well as in-plane magnetic field components. That system of equations characterizes the Alfven-Coriolis 
waves and depends parametrically on the radius: 

- i\v r - 2v e - ] ^ = 0, (44) 

- tXvo + \v r - - \ - ^ = 0, (45) 

2 p z {r)u(r]) dr) 
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_ i\b r - ^ = o, (46) 
drj 

- i\b B - ^ + |S r = 0. (47) 
drj 2 

In addition, the Alfven-Coriolis sub-system (I44p - (I47|I is subject to the vanishing boundary conditions for the in-plane magnetic- 
field components. 

The linear set of equations ()44[) -(147 [ ) ma y be reduced to the following single fourth order ordinary differential equations 
for both b r and bg: 



d_ 

drj 



d j 1 db r e 



+ (3 + 2A 2 )£(r)iL f ) - Axr - 1) n,)/,, „ = 0. 



1/(77) dry 2 \ ^ (77) dry 

Equation (148 1 is the same as the one used by Liverts fc Mondl 1 20091 ) who have derived it for a model problem under the 
assumption of zero radial variations of the perturbations. Here, it is important to emphasi ze however that the radi al coordinate 
is a parameter a fact that renders the radial dependence of the perturbations arbitrary. iLiverts fc Mond (2009) have solved 
(|48p with the aid of the Wentzel-Kramers-Brillouin (WKB) approximation for u(r/) — exp(— rj 2 /2). Remarkably, however, that 
a full analytical solution of Eq. (|48|) is possible for a slightly modified density profile. The first and main step towards that 
goal is to replace the isothermal density vertical steady-state distribution v{rj) = exp(— rj 2 /2) by the following function: 

v(rj) = sech 2 (6r/), (49) 

where the shape parameter b is determined by the requirement that the total mass of the disc does not change, namely: 

f°° 2 f°° 2 

/ exp(—r7 /2)d7j = / sech (brjjdrj. (50) 
Jo Jo 

The result is: 

b = s/2/tt. (51) 

The comparison of true and model number density profiles is presented in Fig. 1 and demonstrates a close correspondence 
between them. Furthermore, it will be shown below that the results derived by employing that profile are weakly distinguishable 
from the WKB results obtained for the true exponential distribution. In fact, notwithstanding the use of the terms true and 
model profiles, such a change of the equilibrium profile of the number density (that is determined by the axial momentum 
bal ance equ a tion) may actually represent some true equilibrium that is obtained from a slightly different gravitational potential 



Spitzer (1943), where a similar density profile has been obtained as an exact solution for flat disc-galaxies whose disc 



mass content is larger than the mass of the central object]. 

Instead of Eq. (|48[) for the perturbed magnetic field it is more convenient to consider the equation for the perturbed 
velocity. Thus, substituting (|49[) - (|51[) into (|44[) - (|47[) . yields the following ordinary differential equations for v r and vg: 

d 2 ( 1 d 2 V r ^\ 2 ^ _ld 2 Vrj) .2, .2 ■,\a2- n 

+ (S + 2\ )p z —- + A (A -±)PgV T ,e=Q. (52) 



i>(rj) drj 2 \v(rj) drj 2 J ^(v) drj 2 

The next step now is to introduce a new independent variable £ = tanh(&J7), such that —1 ^ f ^ 1. Doing that, a simpler 
equation emerges, which may be cast into the following form: 

(L + K~){L + K + )v 6 = 0, (53) 
where L is the Legendre operator of second order: 

L = i [{1 ~ e) i ] > K± = 7^[3+2A 2 ±\/9+16A2]. 

Imposing now the zero boundary conditions of bg at rj — > oo leads to the requirement that the solution of Eq. (|53l) for vg 
diverges polynomially at most when rj — > oo. It is concluded therefore that vg is proportional to the Legendre polynomials 
Pk(£) which are the solutions of the following equation: 

i\ ( (1 ~ c2)£ lr) + K ±Pk ' K± = k{k k = lj 2 ' (54) 

Thus, the eigenvalues A are now determined by the dispersion relation 

J& = %|^[3 + 2A 2 ± v/9 + 16A 2 ] = k(k + 1). (55) 
Employing now the recursive relations for the Legendre polinomyals: 
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\ exp(-i7 a /2) 




Figure 1. The comparison of true, u = exp(— J? 2 /2) (solid) and model, v = sech 2 (br]) (dashed) number density profiles; b = ^/2/tt. 



setting the arbitrary amplitude of ve to unity, and using the linear equations (|44[) - (|47[) result in the following expressions for 
the eigenfunctions that are determined up to an arbitrary radius dependent amplitude factor: 



2 P z (A?) 2 



bt = ^±^[(At-l)(Aj+l)&-3] [ft-i(0-ft+i(0], 
~± fc(fc + 1) ft (1-A*)(A* + !)&-! rn 

o a — — — — — [n-i(U-^+ilg] 



2fc + l 4 



(56) 



where fc = 1, 2, ... plays the role of axial wave number, — for £ = ±1, since -Pjfe(l) = 1 aim ^fc(~l) — (— l) fc f° r au 
Turning back to the dispersion relation (58), it may be written as follows: 



(A ± ) 4 /3f - (A ± ) 2 /? Z (/3 Z + 6) + 9(1 - /§,) = 0, p z 



Pz 



p 



(fc)- 



(57) 



It is thus obvious that the fc-th mode is destabilized when the beta value crosses from bellow the threshold that is given by: 



(58) 



As a result, a universal (for all values of P z (r) and fc) criterion for instability emerges which reads: p z (r) > 1. Written in terms 
of the scaled plasma beta p z the dispersion relation (|57[) has the following solutions for the eigenvalues of the Alfven-Coriolis 
modes (see Fig. 1): 



Pz + 6 ± \J(Pz + 6) 2 - 36(1 -p z ) 



2f3 z 



(59) 



The two eigenvalues, A + and A - , represent fast and slow Alfven-Coriolis waves. While the fast Alfven-Coriolis modes are 
always stable, the number of unstable slow modes signified by the plasma beta. The eigenvalues of the slow Alfven-Coriolis 
modes, A - , are given therefore by: 



A- = a: 



Pz + 6- V(A+ 6 ) 2 -36(1 -AO 
2p z 



, Im{A~) = for p z < 1, 
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+A* ±r~ 




Figure 2. Growth rates ±r a (dashed lines) for the unstable Alfven-Coriolis (MRI) modes and frequencies ±A a (solid lines) for the 
Alfven-Coriolis oscillations vs universal scaled plasma beta, fi z = = Bcr\ 0c? = 2fc(fc + l)/(37r) for the model number density 

v = sech 2 (-^/2/7r?7); k = 1, 2, . . . is the axial wave number. Meshed straight-line asymptotes at j3 z 2> 1 are the scaled Keplerian frequencies, 
±A± = ±1. 



A = iV~ = ±i 



(P z + 6) 2 + 36(/3 z - 1) - z 
2^ 



Re{A") = for p z > 1. 



(60) 



The eigenvalues are imaginary, and the system is spectrally unstable if p z > 1, and real for p z < 1 in which case the system 
is stable. In particular, the minimal critical unsealed plasma beta that is needed for instability is determined by the first 
unstable slow Alfven-Coriolis mode, k = 1, and is given by Per = 0.42. The unstable modes are the well known MRIs and 
from Eq. (|58|l it is easy to calculate how many of them are excited for a given value of the plasma beta. Thus, there are k 
unstable modes for /3cr^ ^ P z {r) ^ fiir + ■ In particular, Eq. ([58} yields approximately the square-root law for the number of 
the unstable modes vs plasma beta: 



k < 



y / 3^W/2. 



(61) 



Relation (|58|) or its simplified version (|61l) is significant for the consequent modeling of non- linear development of the instability. 
It is finally emphasized that the stability criterion as well as the number of unstable modes depend on the radius. Thus, different 
areas within the disc may be characterized by different stability properties as well as different number of unbstable modes. 

The family of the fast Alfven-Coriolis modes, A + (/3 z ), is characterized by the frequencies that are much larger than the 
Keplerian frequency for small values of the scaled plasma beta (large values of axial wave number or small plasma beta) and 
tends to the Keplerian value at large plasma beta values: 



a: 



\ 



Pz + 6+ \J{p z +6) 2 -36(1 -p z ) 
~W Z ' 



ImjA 4 



0. 



(62) 



Expressed in terms of the scaled plasma beta, a single figure depicts all possible stable as well as unstable modes. This is 
shown in Fig. 2. The maximal growth rate for the unstable modes is achieved around p z w 3, which for a given plasma beta 
value determines the axial wave number of the fastest growing modes. 

Frequencies and growth rates for the first five fast and slow Alfven-Coriolis modes are presented in Fig. 3 (in terms of 
the unsealed beta, p z ). The corresponding growth rates F^ obtained by Liverts & Mond (2009) by employing the WKB 
approximation are weakly distingu ishable from tho s e in F ig. 3 and therefore are not shown on the figure. In particular, the 



critical plasma beta Pil found bv lLiverts fc Mondl ([2009) (after correction by factor 2 due to the different definition of the 



plasma beta there) is 0.417, which is practically the same as the one obtained in the current work. 
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Figure 3. Frequencies ±A a for the first five fast and slow Alfven-Coriolis modes and growth rates ±r a for the unstable (MRI) slow 
Alfven-Coriolis modes vs unsealed plasma beta Z for the model number density v = sech 2 {^/2/irri); p£) = 2k(k 4- l)/(37r) is the critical 
plasma beta, k = 1,2, 3, 4, 5 is the axial wave number. 




Figure 4. The toroidal component of the perturbed magnetic field bg vs self-similar axial variable rj = z/H(r) for the Legendre 
polynomials with the axial wave numbers (a) k = 1, 2, 3,4 and (b) k = 6, 15. All curves are calculated for the fixed value of j3 z = 1.5 
Z = Pz/Pcr^ ; Pc? = b 2 k(k + l)/3, dashed and solid curves correspond to the odd and even k, respectively). 



An illustration of the perturbed toroidal magnetic field is presented in Fig. 4. The perturbations are indeed localized 
within the effective height of the disc which weakly depends on the axial wave number. This corresponds to a finite distance 



betwe en the turning points - the natural characteristics of the problem solution in the WKB approximation Liverts fc Mond 



( 2009h ] . To explicitly demonstrate the finite size of the region of perturbation's location, the differential equation for the 
perturbed magnetic field be in the self - similar variable rj is written out: 

^l + 2b^+b 2 K ± (l~£, 2 )b e = 0, f = tanh(&77), 6= A K ± =jfe(jfe + 1), fc = l,2,..., (63) 

which is subject to zero boundary conditions at infinity. Equation (|63[) is derived by using differential equation (1531) for vg 
and a fact that bg up to a constant coefficient of proportionality is equal to dvg jdr\. Transforming as in WKB approximation 
the dependent variable bg — Q(rf) J n(ri)dri, and setting fi(rj) — — £ = — tanh(&77) in order to eliminate terms with first order 
derivatives yield: 



±Q+^( V )Q = 0, ^ 2 ( ?7 ) = b 2 [^ ± (l-^)- l;. (0 1) 



Then the turning points are determined from the relation >c (r\) = 0, and may be found explicitly substituting 1 — £ by 
sech 2 (brj): 

n* = b _1 arcsech — , k = 1, 2, . . . . (65) 

y/Hk + 1) 
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Table 3. Distances of the turning points from the midplane vs axial wave number. 
k 1 2 3 4 5 6 15 35 155 
»7» 1.1 1.9 2.4 2.7 3.0 3.2 4.3 5.3 7.2 




The results for several axial wave numbers are presented in Table 3. Although the oscillation's region unboundedly expanded 
with rising axial wave numbers, this occurs very slowly and the region has relatively low sizes even at high values of k. Note 
that the unbounded growth which corresponds to high wave numbers should be suppressed by dissipative effects neglected in 

the present analysis. 

In order to compare the current results to well known results for infinite homogeneous cylinders iBalbus &i Hawlevl (|199lf ) 
the growth rate of the unstable modes is depicted in Fig. 5 as a function of the axial wave number K, where 

L a k 



K 



HV2 



(66) 



L a = V a /£l is the Alfven length scale, z = Hy2 is the effective height of the diffused disc at which the equilibrium number 
density, n ~ exp[— z 2 / (2H 2 )], falls by factor e _1 . The effective wave number K is the discrete thin-disc analog of the continuous 
wave number for infinite cylindrical discs. For fixed value of the local plasma beta, j3 z = 0.41, 0.5, 1.5, 2.5, 500, the discrete 
set of the points in the plane {K, T^} is presented by one of the interpolating curves 1,2,3,4,5, respectively. The number 
of the discrete points on each interpolating curves corresponds to the admissible values of the Alfven-Coriolis mode number 
k = 1,2,3, ... for which < } < j3 x (r). Also, the range of unstable fc-values is widening as the value of Z is increased. 
In particular, at large plasma beta the correspon ding set of the poin t s due to their large number (curve 5) should tend to the 
continuous curve for infinite cylinder geometry in lBalbus fc Hawlevl (jl99l|). For finite plasma beta there is a discrete number 
of points on each curves in Fig. 5, where the left bound corresponds to the first Alfven-Coriolis mode, k = 1. Thus, for beta 
values close to p„ , the number of unstable modes is small and the disk stability properties significantly deviate from those 
predicted by the infinite cylinder model. 



3.4 Linear stability problem for magnetosonic modes 

For simplicity it is assumed now that the poloidal component of the equilibrium magnetic field is zero and hence j3 z = oo. 
Thus, substituting the exponential anzatz (|42fl into Eqs. (|40|) - (|41|) for the magnetosonic modes, and assuming that their 
frequencies are different from the eigenvalues of inertial-Coriolis waves yield: 
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-i\v z + ±[^-]=0, (67) 
dr) 1/(77) 

1/(77) d*? ^( r ?) dr) z 

Since the effect of the magnetic field is dropped from the linear system (|67[) - (J68J) , the latter actually describes sound waves 
that propagate vertically in the non-homogeneous disc. Eliminating the perturbed number density from (|67p - (|68p in the 
resulting equation yields a second order differential equation for the perturbed axial velocity: 

d 2 v z 1 dv(r\) dv z f 2 d r 1 dv(r))l 



dr] 2 ^(77) dr] dr] I dr) L 1/(77) df] 
or substituting the equilibrium number density, ^(77) = exp(— r) 2 /2), 

lv" 7 *7 + (A - 1)u * = a (70) 

Equation (|70l) is the Hermite equation and the requirement that its solutions diverges polynomially at most when 77 — > ±00 
determines the eigenvalues to be: 

A = ±A S = ±Vm+l, (71) 
while the eigenfunctions are given up to an amplitude factor depending on radius by: 

v z = ^iVrrT+lH m ( V ), i> = e- v2/2 H m+1 { V ). (72) 

Here H m (r)) ( m — 1> 2, ) are the Hermite polynomials. The axial velocities are polynomially unbounded functions at 77 — ^ ±00. 
Note that now the number of the Hermite polynomial m plays the role of the axial wave number for the magnetosonic mode that 
is different from the Alfven-Coriolis modes with the number of the Legendre polynomial as the axial wave number. Formally, 
the above family of the solutions of the eigenvalue problem should be completed by the following degenerate steady-state 
solution: 

v = e~ v2/2 , v z = 0. (73) 
However, the latter may be considered as a modification of the unperturbed equilibrium solution. 



4 DYNAMICAL EQUATIONS FOR PERTURBED THIN DISCS: TYPE II EQUILIBRIA 

Consider now the thin disc problem under the effect of a dominant toroidal component of the equilibrium magnetic fields (see 
Tables 1 and 2). The solutions of su ch problem in the pure hydrodynamic case of zero magnetic field for adiabatic thin discs 
have been obtained and discussed in lshtemler et aL (2010) and are extended here for locally isothermal magnetized discs. 



4.1 The reduced nonlinear equations 

At this stage it is convenient to rescale all the physical variables, /, with the small parameter e, without splitting the solution 
to its equilibrium and perturbed constituents. Thus, in a most general way, the physical variables are given by: 

f(r,(,t) = e S 'f(r,(;,t). (74) 

Here the auxiliary variables denoted by tilde are introduced, and the same orders, S', in e are assumed for the physical 
variables as for the disturbed variables. For the dominant toroidal equilibrium magnetic field such choice is justified by the 
fact that any disturbed dependent variable is of the order of or larger than the corresponding equilibrium variable i.e. S' ^ S 
(see Tables 1 and 2). This imposes an additional limitation on the admissible time of perturbation growth. Substituting (|74[l 
into ([2j - yields up to the terms of the higher order in e the following in-plane components of the momentum equation: 

^r+V r ^JL+V z ^-^ = -^, (75) 
dt dr iK 

dVg ~ dV e T % dVg , , , „ 

- d T +V r-d T -+ V *-dC+ — = > ( 76 ) 

dV z y m y dVz _ c% dh 80 1 d{B 2 e + B 2 r ) m) 
dt dr dC n d( d( 2/3h d{ ( ' 

dh l d{rnV r ) d(nV z ) 

-dt + r^r— + ^T = ' ( ' 







r 


dr ' 


VrVe 


= 0, 


r 
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^1-^=0 (79) 
dBe 8E r dE z _ 

~df + ~dc: dT~ ' (80) 

^ + i^=0. (81) 
at r or 

Here 

E r — V z Bg — VgB z , Eg — V r B z — V z B r , Ez = VgB r — V r B$. (82) 



The system (|75|> - (|82|) is subject the boundary conditions at infinity of least possible divergence for the axial velocity and the 
vanishing of the in-plane magnetic-field components. 



4.2 The linearized problem 

We start by linearizing the MHD equations (|75[) - (|82|) about the steady-state equilibrium solution (|12|l - (|15[) . Substituting 
the decomposition of the total disturbed variables (|18[) into Eqs. (|75|) - (|82|) . yields the following linear set of equations for 
axisymmetric perturbations: 

r)V' 

^ - 2Q(r)V; = 0, (83) 

^ + |fi(r)tf = 0, (84) 
and 

^ + c ^iX ] + lM^!- (85) 
dri d(nV;) 1 c>(rnV;') 

¥ + ~"^^T' (86) 
as;. = ay; 

^r = Bz(r) ^r (87) 



<9t dr 9r 

As for the equilibria type I, the equation for the perturbed magnetic field, B' r , is obtained by differentiating by the equation 
for the magnetic flux function, while the corresponding equation for the perturbed axial magnetic field, B' z , decouples from 
the rest of the equations, and drops out altogether from the governing system (183ll - (|88| l, As in the case of comparable axial and 
toroidal magnetic field components (type I equilibria), also the reduced linear system (|83[l - (|88|l is decoupled into two linear sub- 
systems that describe now the dynamics of the inertial-Coriolis (IC), {V^,v'g}, and magnetosonic (MS), {V' z , n! , B' r , B' e , B' z }, 
modes. A notable difference between the two cases however is that, while for the type I equilibria the radial derivatives drop 
out altogether from the reduced nonlinear system of equations, in the current case, some of the radial derivatives do survive 
the asymptotic procedure, due to the relative smallness of the axial magnetic field component. As will be shown later on, this 
difference turns out to be quite significant as the radial derivatives are responsible for the non modal algebraic growth of the 
initially small perturbations. 

Back to the modes of wave propagation, the IC waves are described by Eqs. (|83|l and (|84p while the MS modes are 
given by the solution of the homogeneous parts of Eqs. (185[) - (|88[ ). Also, the perturbed radial magnetic field, B' r , is completely 
determined by the perturbed radial velocity, V r ', and hence Eq. (|87[) is separated from the rest of the magnetosonic sub-system. 
The right hand sides of Eqs. (|86[) and (|88|l result in the (resonant as well as nonresonant) forcing of the MS waves by the IC 
modes. As will be seen later on, this occurs due to the rotational shear. 

Before actually solving the above linearized system of equations it is noted that, guided by the steady-state solution, it 
is convenient to introduce the following new independent variables: 

9 = *' P= f Jfh> " = ^ = TTTV)' ( 89 ) 

Jo H(r) H(r) fi(r) 

Then derivatives in new and old variables are related as follows: 
d__d_ d__]_d_ jl_JL_J_^ jl 

dt~ 80' d(~Hdv' dr~dp H 2 dp V dr,' [ ' 



Equations (|83l l- (|88[ ) may be rewritten by introducing the following scaled variables: 
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v( /O ,7 7 ,0) = 7 ^t, u( P ,r,,e) = ^- h( P , v ,e) = -^- (91) 
c s (r) N(r) B e (r) 

Below for convenience and with no confusion we revert to the notation t for the new variable 8. This yields the following 
system of equations: 

hw^ = °> <*> 

and 

Q dt dri l v(v) Pev(v) dr] ' 1 ' 

19,1'. dv z dv r = = ,dv r . 

+ W " ^ z = ~W ~ + n( W " ^ H (95) 

Q dt 8(r) dr] ' 1 ' 

1 <%e 9« z dv r 3 1 3ne „ n 

n~a7 + ~cV ~~dp~ ~2 r + It~P)~d^~ ° BVr + D ^ (97) 

Here Pe{p), Pz(p) and S(r) are given by (|35|l ; j3$(p) and j3 z (p) are proportional to the characteristic lasma beta, /3, that is 
now based on the dimensional toroidal component of the equilibrium magnetic field as the characteristic scale, B* = Bg(r*)\ 
Dn{p) , Dn(p) and Db(p) are the following logarithmic derivatives: 

= dln(c s rN) - d]n(c s /n) _ d\nH - d\n{c s B e ) 

L> N = , Dn = = — — , D B = j ■ (98) 

dp dp dp dp 

A single equation for the magnetosonic modes may be derived now by eliminating the toroidal component of the perturbed 
magnetic field as well as the perturbed number density from Eqs. (|94[) - (|97[1 . The result is: 

1 d 2 v z . 1 .d 2 v z dv z , , , ,..„., ■ , ■ , . 

i2 2 dt' 1 Pevyn) orf dr] Pevyn) 2 or] ]j /3 Z dr] 1 fievyn) dr]dp 

where L z is the linear ordinary differential operator: 

= dv r = r , 1 9 9u r d(r] 2 v r ) - 1 dv r 

L z v r = Djv- D n (1 + - )■— ( 7? — -) !_ — i + Da — . (100) 

dr] Pev{r]) dr] dr] dr/ Pe^in) dr l 

Note that in the important particular case of pure hydrodynamic discs the equilibrium and perturbed components of the 
magnetic field should be set to zero, i.e. f3g = j3 z = oo and b r = bg — 0, then Eqs. (|97|l - (|98|) with the properly resolved 
uncertainty in Pe/j3 z are satisfied identically, and Eq. (|99|l is replaced by 

1 d 2 v z d 2 v z dv z d 2 v r = dv r = 9 d(rfv r ) 

+ 1- v z = » » + Div — ^M^W - ^ - ) Ha — ]■ ( 101 ) 




f2 2 cfe 2 c^ 2 97; 9r;9/9 dr] dr] dr] dr] 

The above problems formulated for the magnetized and magnetic-field free discs are complemented by the boundary 
condition of least possible divergence at infinity for the axial velocity, and in the case of magnetized discs the vanishing 
conditions for the in-plane magnetic-field components should be also satisfied. 

An important result is immediately seen from the above set of linearized equations: to leading order in e: Eqs. (|92|I - (|93H 
decouple from the rest of the system in the thin disc approximation. Note additionally that Eq. (|96l) for the radial magnetic 
field, explicitly determined by the radial velocity, may be solved separately from the rest Alfven-mode sub-system. Thus, the 
sub-system (91)-(92) describes pure hydrodynamic in-plane inertia- Coriolis waves, and signifies the decoupling of the latter 
from the Alfven waves. This occurs due to the negligibly small values of the projections on the disc plane of the pressure 
gradient and Lorentz force in the momentum balance equations (of the order of e 2 ) compared with the inertial terms (of the 
order of e°) in the thin disc approximation. Thus, in the case of poloidal-dominated equilibrium, the magnetic field is of order e° 
which renders the Lorentz force of the same order of magnitude as the inertial terms, and thus couples the Lorentz and inertia 
forces. This results in the coupling of the Coriolis and Alfven waves that eventually leads to the MRI. In contrast, the case of 
the toroidal-dominated magnetic field inevitably results in the scaling B z ~ e. This, as indicated above, means that the radial 
and azimuthal components of the Lorentz force are negligible with respect to the corresponding components of the inertial 
forces in the radial and azimuthal components of the momentum equation. As is already apparent, that decoupling between 
the inertial and Lorentz forces causes to a crucial deviation from the former case (i.e., with polidal-dominated equilibrium 
magnetic field) and leads to the removal of the mechanism that is responsible to the MRI. 
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4.3 Linear stability problem. 

The perturbations may be presented as follows 

f(p, n, t) = exp[-iAO(p)i]/(p, 77), (102) 

where A = ui/Q. is the scaled frequency. 

Inertia-Coriolis modes. Substituting (1 102 p into (|92|l - (|93[) yields: 

- i\v r - 2v = 0, (103) 

- i\v e + \vr = 0. (104) 
The corresponding dispersion relation is therefore given by 

A = ±l, (105) 

which represents, the stable epicyclical oscillations in the disc plane, since for Keplerian rotation, the epicyclical frequency 
equals to Keplerian one, x(p) = &(p)- As p and n are mere parameters each ring p = const, n — const vibrates independently 
in its own plane. If viscous stresses are taken into account, an axial profile is imposed due to the mutual shear stresses 
between the rings and each entire cyli n drica l s hell vibrates independ ently. This is indeed the case that has been solved in 
Umurhan et~aiT(|2006h ; iRebusco etahl l|200gh ): Ishtemler et all |201oh ]. As the axial and radial coordinates play the role of 



passive parameters, the eigenfunctions of the inertia-Coriolis modes are determined up to arbitrary amplitude A(p, n) from 
<P3>(fTU5l) : 

ve = TiTjVr = ^fi-A(p,n). (106) 

Since any special form of the function represents some given set of initial conditions, the following self-similar separable form 
of the planar velocities is considered as an example: 

A(p, V ) = F(p)G( V ), (107) 

where the F(p) and G(n) are arbitrary functions. 

Magnetosonic modes. For frequencies that are different from the eigenvalues of inertial waves, namely A = cu/il(p) = ±1 
in (|105[) . both the perturbed in- plane components of the velocity as well as the radial component of the perturbed magnetic 
field are zero, i.e. v r = i>e = b r = 0. Therefore, the perturbed axial velocity, number density, as well as the toroidal magnetic 
field are described by the set of equations that governs the dynamics of the magnetosonic waves: 

/3 g v(rj) dn 1 dn 

- lA -^T + ^-Vv z = 0, (109) 
u(n) dn 

-i\b e + ^ = 0, (110) 
dn 

is subject to the boundary conditions at infinity of the least possible divergence for the axial velocity and the vanishing 
conditions for the in-plane magnetic-field components. Partial solutions are presented below . 

(i) The pure hydrodynamic system. First start with a simple exact solution for the pure hydrodynamic system. Accounting 
to Eq. (fTOTjl this yields 

rj— h (A - l)v z = 0, (111) 



dn 2 dn 

-i\J^ + ^-r)v z = Q, (112) 
u(p) dn 

complemented by the boundary condition of least possible divergence at infinity for the axial velocity. That problem is 
degenerated to the sound waves, and coincides with the case of the poloidal-dominated equilibrium magnetic field, and its 
solution is: 



A = ±A S = ±Vn7+l, v z = -i\W(p)H m (ri), V - = W{p)H m+ i{v) at 77 = ±00. (113) 

Here A are the radius-independent eigenvalues; W(p) are the radius-dependent amplitude factors; H m (n), m = 0, 1,2, ... are 
the Hermite polynomials. 

(ii) The limit of large plasma beta in magnetized discs. The WKB approximation is employed now. In general, it may 
be applied for a wide range of plasma beta values. However, in order to simplify calculations, we restrict ourselves here by 
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its application to the limit of large but finite values of the plasma beta. Such approach is necessary as a simple asymptotic 
expansion in large plasma beta ft$ fails at large 77 (small v(r))) due to the non-uniform limit of j3$v(rj) in Eq. (|108[) . As a first 
step the following equation for bg is derived from the system (|108[) - (|110[) 

d% Wev(n) - 1 dbg 2 p e v{ri) t , 

W ~ *^(„) + fc*r + (A " TT^Fw = °' (114) 

with the following boundary conditions: 

b = at 7/ = ±00. (115) 

To reduce Eq. (|114[) to the form appropriate for WKB approximation, the dependent variable is transformed in the following 
way in order to eliminate terms with first order derivatives: 

b e = Q f^ v , ^-i- fy;:; . (ii 6) 



2Pe Pev{n) 



Then 



0+x 2 (r ? )Q = O ; (117) 
where 



*(„) = - 2 1 + ^ gagW _ 1 _ 2 1 l-^ 2 ^) 2 a _ ^£(g) 

W ^ 2/3 e [l + ^r,)] 2 Zfal+pePin) 4/3% 1 + P~ev[rf) K 'l + /3 e u(v)' 

The approximate solu tions appropri ate to real and imagine >c(r\), at which the solution oscillates and decrease, respectively, 
are as follows (see e.g. iMigdall (|1977l l): 

W f W f 

Q = — exp(+7. / H\7\)dr\ and Q — — exp(+ / | x[rf) \ drj. (119) 



The regions of each of the solutions given in Eq. (|119[) are determined by the turning points which are the zeros of x{rj). It can 
be easily seen that large-plasma beta approximation is uniformly valid in the region between the turning points, and violated 
at large 77 due to uncertainty in the term f3$v(r)) in (|118[l for exponentially vanishing v{rj) = exp(— rj 2 /2) at 77 — > ±00 and 
Pe ~ P 2> 1. Since application of the WKB method is justified asymptotically at large frequency A, Eq. (|118|) is considered 
for large but fixed values of j3g and A, and for such large 77 that both j3gi/(rf) and j3gp(r)) are small. Then Eq. (|118[) yields to 
leading order 

S(v)m-±-(1 + £-), KV) (120) 
2pg Ipe Ape 

Choosing the minus sign to satisfy the boundary condition (|115[) in the approximate solution (|119[) for 77 |>| 77, |, yields the 
following form for sufficiently large 77: 

i 1 , V 2 , f n 1 1 . v We , if 

be ~ -, r oxp( =- ± / x 1377) = — — exp( = 

M Py 4/3 fl A, 1 1 " 77 2 PV 2f3g 



■cxp(-^±/ |x|d77) = ^exp(-^). (121) 



To further simplify the calculations, assume additionally that fie S> 1, then the large beta limit is uniformly valid within the 
interval | 77 |^| 77, |, and >c and \i are given by 



2/. \ 1/2 2s / n 1 



K\ri)n-~(r, *-r£), ^^-77, (122) 



where 77, = ±\/4A 2 — 2 are the turning points of the equation (|117|) which separate the regions of oscillatory and exponentially 
decreasing behavior of the solution. Thus, h = at | 77 | = | 77, |, x 2 is positive within the inner regionj 77 |^| 77, | and negative 
otherwise within the outer region. To completely determine the problem solution, the conventional WKB approximation applies 
the matching conditions to the solutions (| 11 9p in the inner (| 77 |^| 77, |) and outer (| 77 |^| 77, |) regions in order to obtain the 
coeffic ients in (1119)1 and the eigenvalue equation. The result is the following Bohr-Zommerfeld condition Landau fc Lifshitsl 
|l997l) ] that determines the dispersion relation 

+v, 1 

xdr] = 7r(m+ -), (123) 
-v* 2 

where x is given by (|122|) . and m is the number of zeros of the solution within the inner region. Substituting (|122|) into (|123[) 
yields 

A = ±Vm+T (124) 
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The eigenvalues (|124[) exactly coincide with those in (|113[1 obtained with no restrictions on the value of axial wave numbers 
m for the pure hydrodynamic system. 

Finally note that although the WKB method is valid asymptotically for large frequencies, A > 1 (i.e. large plasma beta 
or, alte r nativ ely, large axial wave numbers, m), it frequently yields a fair approximation up to finite values A ~ 1 [see e.g. 
Migdall (1977)]. This allows expecting that eigenvalues (|124|l will be appropriate even for low axial wave numbers m = 1, 2, . . .. 



(Hi) The limit of small plasma beta in magnetized discs. Consider the problem (|108[) - (|110|l in the limit of small plasma 
beta /8e(p) ~ /? <C 1. Assuming simultaneously that the eigenvalue A is large, so that \ 2 f3g ~ f3° , Eq. (|108[) yields to leading 
order in small ft 

+ \ 2 l3oHv)vz = 0. (125) 

Replacing as in Section 3.3 the steady-state isothermal density distribution v(r]) = exp(— if /2) by the hyperbolic function 
i>(ri) = sech 2 (b?7) with b = and introducing the auxiliary variable £ = tanh(&7/), transform Eq. (|125[) to the following 

form: 

Lv z + \ 2 v z = 0, (126) 
where L is the Legendre differential operator of second order, 

and A 2 is given by: 

A 2 = ^A 2 (p)A,(p)~^°. (127) 

Imposing now that v z diverges polynomially at most when <r\ — ¥ oo, leads to the conclusion that the solution of Eq. (|127[) for 
v z is proportional to the Legendre polynomials P m (£): 



P m (0~P , A = m(m + 1) ~ p, m = l,2,.... (128) 



W{p) 

Equations (|106p - (|l(J7p yield that bg and v are of higher order in small plasma beta: 

?V 2 ~, fl!/2 



W(p) 



(1 - e)[^ + ^-arctanK)P m (C) - P m -i(t)W ~ ^ ■ ( 12 9) 



W(p) b]j m + l K mb 

Since P m (l) = 1 and P m (— 1) = (~l) m f° r a U m , bo satisfies the vanishing boundary condition at infinity r\ — > ±oo (£ — > ±1). 

Thus, summarizing, the main conclusion from both limits (high and low plasma beta) is that the inertia-Coriolis waves 
are decoupled from the Alfven waves, which henceforth leads to the complete stabilization of the MRI's. This, as will be 
subsequently seen, leaves the stage exclusively to the non modal algebraic growth mechanism. Remarkably that both limits 
assume sufficiently large values of the eigenvalues (scaled frequensies A): in the limit of large /3 in (|124p due to high frequency 
assumption inherent to WKB, while in opposite limit of small j3 in (|128p . because to asymptotic behavior of A ~ /3 _1//2 . 



4.4 Non-resonantly and resonantly driven magnetosonic modes 



There are two mechanisms responsible for algebraic time-growth of perturbations, namely non-resonant and resonant excitation 
of vertical ma gnetosonic wave s by p l anar inertia-Co r iolis m o des. This fact for the pure hydrodynamic adiabatic discs was 
established in [Umurhan etafl (|2006h : iRebusco et all jiooil )): Ishtemler et all |2010[ )]. As shown below, both mechanisms are 



also relevant in a modified form for vertically-isothermal magnetized discs. Due to the rotation shear effect and the presence 
of radial derivatives in the system (|95p - (| 10 1 j) . the amplitude of the magnetosonic modes driven by the inertia-Coriolis modes 
grows linearly in time in the case when the magnetosonic mode frequency is not equals to the inertial-Coriolis one. While for 
the same frequencies of the magnetosonic and inertial-Coriolis modes, the latter grows quadratically in time. 

To demonstrate linear and quadratic in time growth of the perturbations, both non-resonantly and resonantly driven 
modes are illustrated below by simple explicit solutions. The non-resonantly driven magnetosonic waves are considered in the 
limit of small toroidal plasma beta and for pure hydrodynamic system, as typical examples. Noting that the resonantly driven 
magnetosonic modes may exist for arbitrary frequencies, we restrict ourselves by considering a pure hydrodynamic system 
which is characterized by simple explicit solutions with no limitation on the frequency value. 

The driving modes are characterized by Eqs. (|102p and (|106p . As a result, the driven magnetosonic or sound modes in 
magnetized or magnetic-field free discs, i.e. (|99p or (|101l) . respectively, are described by the following expressions: 
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{v r , v$, b } = exp[-i\Cl(p)t]{v r , v e , b e }(p,rj,t). (130) 

Inserting then Eq. (|130[) into the system of equations ( (|94l) -f (| 1011) and keeping only terms that are proportional to t on their 
right hand sides yield 



H 2 v z 1 * 2 



&-w t 1+ ^^ + ^ +0 - =Tft ^ [ -^ wtI(1+ ]dw ) ^W' (131) 

kl [ lk) ] + ^~ VV ^ ^^(PW f/r, (132) 

1 db e dv z . =. 

n aT + ftT = ±**<*p[t*W]^. (133) 

The problems ()131I) - (I133|) is subject to the boundary condition of least possible divergence for the axial velocity, and for the 
magnetized discs complemented by the vanishing conditions of the toroidal magnetic field at r\ — ±oo. As is indeed seen from 
that system, the magnetosonic wave is driven by the inertial wave, while the linear growth in time is entirely due to effect of 
the rotation shear, and disappears for dCl/dp = 0. 

As was mentioned above the two limits of small toroidal plasma beta as well as pure hydrodynamic system are considered 
analytically. If the scaled eigen-frequency of the inertia-coriolis waves (|109[) , A = ±1, do not coincide with the eigen-frequencies 
of the magnetosonic waves (|118fl . A = ±ym + 1, i.e. at m / 0, non-resonantly driven magnetosonic modes are excited. The 
stable magnetosonic modes may be excited resonantly by the inertial modes if any pair of respective eigen- values (|109|) with 



A = ±1 and (| 1 1 8|) with A = ±vm+T, coincide. It is easy to see that this may happen only for m = 0. 

(i) Non-resonantly driven magnetosonic modes in the limit of small plasma beta. At vanishing plasma beta, /3, in the 
leading order approximation Eq. (|131|) leads to a degenerate quasi-steady problem for v z on the Keplerian time scale. Twice 
integrating Eq. ()131[) in r\ and setting arbitrary integration constants to zero, yield: 



v z — ±it-r— / irdrj. (134) 
dp Jo 

The rest unknown number density and toroidal magnetic field can be found at known v z from non-degenerate unsteady 
relations (THE?)) -|fI5g)l. 

(ii) Non-resonantly and resonantly driven magnetosonic modes in pure hydrodynamic discs. Setting /3g = oo, v r on 
the right-hand sides of (|131l) - (|133|) may be expanded in terms of the complete set of the spatial eigen-functions for the 
magnetosonic mode, namely, in the Hermite polynomials: 

MP,V) = W r , m {p)H m (ri). (135) 

m=0 

Examining the effect of a single term in the expansion above and utilizing the following recursion relation for the Hermite 
polynomials: 

dH m _ m jj m ' \ m = l,2, 
at] 

substituting the exponential anzatz (|130[) into Eq. (|101|) . and keeping only terms that are proportional to t on its right hand 
side yield 

1 9 2 «z,m-i d 2 v z , m -i , dv z , m -i .= ,,,(in„, , , rr , . . s 

Q2 Qtf q^2 ^ — q^j r-w*,m-i = ^-it exp[=piS2(p)tJ —W r ,m-i (p)mH m -i(r]) , (136) 

where for m = 1 the right hand side of Eq. (1 136 p is a solution of the homogeneous part of the equation and hence represents a 
resonant driving force. For any other value of m the magnetosonic waves are driven non-resonantly. The typical non-resonant 
( m = 2) and resonant (m = 1 ) solution of Eq. (I136[) may be written in the following form: 

v z ,m-i = Vz, m -i(p, r), t) exp[^fiCl(p)9] = ^fiexp[^fiCl(p)t]W z ,m-i(p,t)H m -i{r]), (137) 
W zA (P,t) = (ttn)W^(p) + W^(p) for m = 2, (138) 

W z Ap,t)^^^W^(p) + (im)W^ \p) for m = l. (139) 

Here the coefficients are determined through an arbitrary amplitude of the radial velocity W r ,i(p) for non-resonant (m = 2) 
and resonant (m = 1) system, respectively, are as follows 

W^(p) .2dQ) W<%) 2dQ 

WApj = - l 7il[p-> WAp)= ±l Tilp- for m = 2 ' (140) 

W<$(p) .ldti) W™(P) .lldti 

-rr-r- — ; — r- = ^2— — , — —. ; = —1- = —— IOT 171=1. 141) 

W r ,x{p) Qdp' W r ,i(p) 2Q. dp y ' 
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The corresponding expressions for number density and toroidal magnetic field may be derived from the rest equations of 
the system Q5g]| -JT33J. 



5 SUMMARY AND DISCUSSION 



A comprehensive asymptotic analysis in small aspect ratio of the disc, e, has been carried out of the response of thin vertically- 
isothermal Keplerian discs to small magnetohydrodynamic perturbations. Two regimes of axisymmetric instability have been 
identified depending on the type of equilibria in small aspect ratio of the disc, e. The first is developed in the axially dominated 
equilibrium magnetic configurations, B z ~ e°, and is excited at a relatively low level of hydrodynamic perturbations in the 
disc plane V£, Vg ~ e 1 , while the second regime occurs in the toroidally dominated magnetic fields, Be ~ e°, at relatively high 
level of hydrodynamic perturbations, Vr, Vg ~ e° (the axial velocity perturbations are in both regimes of the same orders, 
Vg ~ e , see Tables 1 and 2). As a result of that scaling, the pure hydrodynamic limit is achievable only within the second 
regime. Indeed, it is demonstrated that as distinct from the first regime, the second regime can not produce spectrally unstable 
modes, but excites a weak algebraic non-modal growth in time. The inertia-Coriolis driven magnetosonic mode leads to their 
non-resonant and resonant coupling that induces, respectively, the linear and quadratic in time growth of perturbations. 

It should be emphasized though that that result is restricted to axisymmetric perturbations and nonaxisymmetric ones 
may give rise to spectral instabilities in the domi nant toroidal field case a s well as indeed is the case for magnetized Taylor- 
Couette flows (see Rudiger et al. 2007). Indeed, Terauem fc Panaloizou ( 19961 ) have derived sufficient conditions for such 
instabilities in some sense of a thin disk limit. If such instabilities i ndeed exist they are of a t ransient nature as the radial wave 
number grows and drives the system out of the instability regime, iBalbus fc Hawlevl 1119911 ) . In addition, as demonstrated in 
Sections 3 and 4 spectrally unstable perturbations may occur if the Lorentz force is of the same order as the inertia terms. 
This means that the spectrally unstable non-axisymmetric perturbations are expected to be of higher order in e than those 
considered in Section 4, that grow linearly and quadratically with time. Thus, it takes full nonlinear numerical calculations 
to determine which perturbations are more significant, namely, non-axisymmetric ones that grow exponentially from a very 
low level with a decreasing growth rate, or linearly and quadratically growing symmetric perturbations that grow from initial 
disturbances of much higher amplitudes. 

For perturbations of type I the main accent is made on the dominant poloidal magnetic field, while the general case of 
comparable poloidal and toroidal components of the equilibrium magnetic field is discussed in Appendix A. Explicit solutions 
of the stability problem for the dominant equilibrium poloidal magnetic field of type I are obtained by replacing the true 
isothermal density vertical steady-state distribution P(rj) = exp(— rj 2 /2) by the hyperbolic function v(ji) = sech 2 (&r;) which has 
the similar shape and the same total mass of the disc. This model profile represents some true equilibrium that is obtained 
from a slightly different gravitational potential. The model eige nfunctions are e x act so lutions of the model very close to 
those obtained approximately by WKB from the true problem Liverts fc Mond ( 20091 )]. but they form the full family of 
explicit, simple and orthogonal eigenfunctions of the model problem. Such properties are significant for the consequent study 
of non-linear development of the instability. 



A study of the stability of isothermal disks has been presented before in iGammie fc Balbus I (|1994T ) and lLatter et al 



( 2O10h . who have also utilized the model density profile sech 2 ??. The main differences between those two works and the current 

investi gation are: 

(Ti bammie fc Balbus I \ 19941 ) and lLatter et all |201ol ) have approached the stability problem of the disk within the local 
setting of a shearing box described by local Cartesian coordinates. In contrast, in the current work the entire disk is considered 
and is described by the natural cylindrical coordinates with the aid of asymptotic analysis in the natural small parameter of 
the problem, namely, the ratio of the thickness of the disk to its characteristic radius ~ e. 



(ii) Ste mming from t he she aring box local approach, the solutions for the perturbations that are given in lGammie fc Balbus 

! 19941 ) and lLatter et al.l (|201oh are independent of the radial-like coordinate. In the current work, however, the radial deriva- 
tives drop out of the equations due to the asymptotic analysis in small e. Consequently, the radial coordinate is a mere 
parameter in the subsequent calculations and conclusions. Thus, the result is that perturbations are excited on each separate 
cylindrical shell, independently of the other shells, with amplitude that depends arbitrarily on the radius through some initial 
conditions. As a result, it comes out rigorously that the stability conditions depend on the radius, and different rings within 
the disk have different stability properties, like number of unstable modes as well as growth rates. 

(iii) The current work includes also analysis of the acoustic spectrum. Though stable in the case of purely axial magnetic 
field (explicit expressions are derived for the eigenvalues as well as for the eigenfunctions), as is sketched in Appendix A for 
large plasma beta, a toroidal component of the same order of magnitude as the axial one leads to the driving of acoustic 
modes by unstable MRI's. 

(iv) Figure 1 demonstrates that best choice of the parameter b is obtained from the requirement of identical total mass 
of the disc for the two density profiles, v(ri) = exp(— rj 2 /2) and u(ri) = sech 2 (brj). 

In addition a qualitative analysis of the influence of the equilibrium toroidal magnetic field of type I on the disc stability 
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is carried out in Appendix A in the limit of large plasma beta on two characteristic scales of AC's and MS's modes. In that 
case as distinct from the small plasma beta (/? < 1), the AC and MS modes are decoupled. On the characteristic scale of 
AC's mode, this leaves the resulting dispersion relation the same as in the case of the dominant poloidal magnetic field (zero 
toroidal magnetic field), and the influence of the toroidal magnetic field is reduced to excitation of the AC-driven MS mode. 
On the characteristic scale of MS's mode, a stable MS mode decouples from the AC mode and the influence of the toroidal 

magnetic field is reduced to excitation of the stable MS-driven AC mode. 

In the present work the model of the non- modal instability Shtemler et al. (2010)] for pure hydrodynamic adiabatic discs 
is adopted for more realistic vertically-isothermal discs with diffused density vanishing at infinity. The model is also extended 
on magnetized discs with type II equilibria. Within that model the eigenvalue problem has been solved in the characteristic 
cases of small and large plasma beta as well as for pure hydrodynamic systems. Furthermore, the magnetosonic waves driven 
non-resonantly with inertia- Coriolis ones are considered in the limit of small toroidal plasma beta, as a typical example. Since in 
the limits of large and small plasma beta, the magnetosonic modes are excited with high scaled frequency, A > 1, they can not 
be resonantly driven by the low- frequency inertial-Coriolis mode with Ail. By this reason the resonantly driven magnetosonic 
modes are calculated for only the pure hydrodynamic system, in order to illustrate a quadratic in time growth of perturbations. 
The rotation shear is found to be responsible for the linear and quadratic growth in time. Furthermore, compressible inertia- 
coriolis stable oscillations continuously pump energy from the sheared steady state equilibrium and transfer it (resonantly as 
well as non resonantly) to the continuously algebraically growing magnetosonic waves. Compressibility is indeed inevitable 
due to the supersonic rotation of the steady-state disc combined with its small vertical dimensions which make the vertical 
magnetosonic crossing time of the order of a rotation period. As distinct from MRI the non modal growth is exhibited for all 
admissible parameters of the system. 

Note that the present study has been restricted by consideration of the axisymmetric perturbations only. In particular, 
for type II perturbations the adopted analysis is considered as the necessary step in achieving a complete picture of transition 
to turbulence in accretion discs. First, it provides the instability for the case of very weak magnetic fields close to pure 
hydrodynamic perturbations in the discs. Furthermore, since the magnetized discs are stable with respect to axisymmetric 
normal modes, an alternative scenario of an algebraic instability has been investigated. The present analysis of both resonantly 
and non-reso nantly driven magnetosonic mo des is interesting in the context of the well known scenario for sub-critical transition 
to turbulence Schmidt fc Henningsonl ( 2001 ) , since the non-spectral mechanism of the power growth in time provides instability 
even for the system parameters of the discs stable with respect to normal modes. The considered scenario leaves to be vali d 
for sub-critical parameters of t he discs within more gener al model of non-axisymmetric MRI [e.g. iBalbus fc Hawlevl (|1991h . 



Terquem fc Papaloizou I |l996l ). iFoglizzo fc Taggerl (|l995l )]. 



APPENDIX A: STABILITY FOR COMPARABLE IN-PLANE EQUILIBRIUM MAGNETIC FIELDS 



The general relations for comparable in-plane components of the equilibrium magnetic field of type I are presented below. 
Substituting the exponential anzatz (|42[1 - (|43[ ) into the general system (|36[) - (|41|1 results in the following system of linear ordinary 
differential equations (parametrically depending on the radial variable r) for the perturbed toroidal magnetic field and axial 
velocity, bg and v z : 



1 



L AC b e = §S(r)^(r){J- 



1 d 2 v z 



1/(77) dr\ 2 



.3(A*-l)tof} 



di] 



S(r) d 



be 



Here S(r) is the coupling coefficient between AC and MS modes, 



LAcbe 



r - _ ^fi- 

LmSVz = —TT + 



d_ r_L <T_ 

dr] I v(rf) drf 2 

1 



1 dbg 
V(ri) dr) 

dv(j]) dv z 



} + (3 + 2A 2 ),3 z (r) 



+ 



d_ 

dr) 

dv{rj) 



]} 



1/(77) dr) 

d 2 v 



+ A 2 (A 2 -l)/3 2 (r)6 e , 



dv z ,2 



1)t5z 



(Al) 
(A2) 

(A3) 
(A4) 



drf ' 1/(77) dr] d-q ' I ' drjlv(ri) dr] J J "* d?7 2 

Equations (|A1[) - (|A2[) are complemented by the boundary conditions of the toroidal magnetic field vanishing and of the least 
possible divergence for the axial velocity at infinity. The operators Lac and Lms are such that the two principal, AC and 
MS, modes satisfy the following homogeneous problems for zero equilibrium toroidal magnetic field (i.e. for zero coupling 
coefficient S(r) = 0) with the same boundary conditions: 



LAcbg - 
LmsVz 



0, 
= 0. 



(A5) 
(A6) 
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Below that problem will be considered for large plasma beta on two characteristic scales of AC's and MS's modes. 
Al Zero toroidal magnetic field. 

First consider the problem for zero coupling coefficient, S(r). Note that if set S(r) = 0, the equilibrium toroidal magnetic field 
has no input in the thin disc approximation under consideration, the problems for AC and MS modes are decoupled, and can 
be considered separately (as in Section 3 above) . The method of two-scale asymptotic expansions in large plasma-beta will be 
applied below, which is convenient for the qualitative analysis of the disc stability. The method is equivalent to WKB approach 
to the problem (IA5I) for MRI mode that also implicitly assumes ft >> f. In that connection note that as is demonstrated 
in Liverts fc Mondl (|2009h the WKB approximation for MRI mode well describes the exact numerical solution up to plasma 



beta /3 ~ 1 (see also Section 3 where the WKB approximation is mentioned very close to the exact explicit solution for model 

density profile). 

According to the method of two-scale asymptotic expansions [see e.g. Navfeh 1 19731 )]. let us introduce the slow, f\ , and 
fast, i), variables at /3 » 1 for AC mode: 



fj = r), v = P a £ fffaK (a>0), 3,(r) = M^~^°, ± = F[g(f})±+p-*±], (A7) 

b e (n) = b e ,o(fj,v)+l3~ a be,i(v,n) + ---, A = Ao + P~ a M + — (A8) 
Applying the principle of the least possible degeneration of the problem yields: 

a =1/2. (A9) 
To make the coefficients of the leading order equation (|A1|) independent of the slow variable, set 

g(v) = v 1/2 (v)- (A10) 
This yields 

+ (3 + 2\ 2 )/3 z (r)^jf- + A|(Ag - l)#(r)8», = 0. (All) 

The vanishing at infinity solution of the linear equation (|A11|) can be easily found explicitly to leading order in ft up to 
amplitude factors that depend on slow variable, fj. The equation (|A11|) evidently conserves all the principle derivatives by fast 
variable in the exact Eq. (|A5[) . The dependence of the amplitude factors on the slow variable are determined by the problem 
of the next order approximation in large plasma beta. To leading order in j3, these two first leading order problems completed 
by the corresponding boundary conditions fully determine the approximate solution of the problem including the dispersion 
relation for the eigenvalue Ao- 

A2 Comparable toroidal and poloidal magnetic fields at the characteristic scale of AC mode. 

Let us now apply the two-scale asymptotic expansions (|A7|) - (|A8[) in ft » 1 to complete relations (|A1|) - (|A4I) for the mixed 
toroidal and poloidal magnetic fields on the characteristic scale of AC's mode. It is assumed as above in the case of separated 
AC mode that a > 0, and additionally 



S W = \¥f\ ~ Mv) = P c m[tizAv,v)+r a t*Av,v) + ■■■}, A = A () + /3- 1 A 1 + .... (A12) 
V Pe( r ) 

(i) Leading order approximation. The principle of the least possible degeneration of the problem yields: 

a = 1/2, c=-2. (A13) 

To make the coefficients of the resulting relations to leading order independent of the slow variable, set 

g(v) = v 1/2 (rj), f(v) = "~ 3/2 (v)- (A14) 
This leaves Eq. (fTTj) valid for the AC mode 

+ (3 + 2\t)/3 z (r)^jf- + \l(\ 2 - l)p 2 z (r)b efi = 0, (A15) 
while the MS mode is the AC-driven mode 

dv z ,o ., S(r) t ... 
= " !A »T7 tVo- (A16) 
9v /3 z (r) 
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Thus, in the limit of large plasma beta, the AC mode decouples on its characteristic scale from the MS mode, and the influence 
of the toroidal magnetic field leaves the problem for the AC mode the same as in the case of zero equilibrium toroidal field, 
and simultaneously induces the AC-driven MS mode. 

(ii) On the next order approximations. Note that the leading-order governing relation for 6e,o is of the order of p 2 , while 
the next-order relation for bo,i is of the order of /3 3//2 . Estimating the terms, which are proportional to the coupling coefficient 
S(r) and induced by the equilibrium toroidal magnetic field, note that they are of the order of /3 1//2 , i.e. much less than f3 3 ^ 2 . 
Hence the toroidal magnetic field has no input neither to the leading- nor to the next-order approximations, which determine 
the resulting approximate solution of the problem. In particular, this leaves the resulting dispersion relation the same as in 
the case of zero equilibrium toroidal magnetic field considered in Section 3. 

A3 Comparable toroidal and poloidal magnetic fields at the characteristic scale of MS mode. 

In that case the conventional asymptotical method of subsequent interactions may be applyied: 



s ^ = \¥f\~P°> &M = ^W, (M7) 

V Pe( r ) P 

be(v) = bg,o(v) + P^beAv) + ■■■■, Mv) = v z ,o(v) + P'^zAv) + •■-, A = A + /3 _1 Ai + . . . . (A18) 
Substituting (A17)-(A18) in (Al) - (A4), and applying the principle of the least possible degeneration of the problem yield that 



the equation for v z ,o and &e,o( r ?) are reduced to the conventional relation (|A6f) for MS mode with the eigenvalues Ao = V rn + 1 , 
(m = 1,2,.. ., see Eq. ((H) ) 

n— h (A - l)Vx,a = 0, (A19) 



.o = -^(r)=^. (A20) 



drj 2 ' dn 
and to the relation for the MS-driven AC mode 

2\f {r) d V 

Thus, on the characteristic scale of MS's mode, the stable MS mode decouples from the AC mode at large plasma beta, and 
the influence of the toroidal magnetic field is reduced to excitation of the stable MS-driven AC mode. 

Finally note that above qualitative asymptotic analysis of the problem in the limit of large plasma beta demonstrates 
both spectral and transient stability of the disc embedded in the mixed equilibrium toroidal and poloidal magnetic fields of 
comparable values. 
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